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We consider several spin-unrestricted random-phase approximation (RPA) variants for calculating correlation energies, 
with and without range separation, and test them on datasets of atomization energies and reaction barrier heights. We 
show that range separation greatly improves the accuracy of all RPA variants for these properties. Moreover, we show 
that a RPA variant with exchange, hereafter referred to as RPAx-S02, first proposed by Szabo and Ostlund [A. Szabo 
and N. S. Ostlund, J. Chem. Phys. 67, 4351 (1977)] in a spin-restricted closed-shell formalism, and extended here 
to a spin-unrestricted formalism, provides on average the most accurate range-separated RPA variant for atomization 
energies and reaction barrier heights. Since this range-separated RPAx-S02 method had already been shown to be 
among the most accurate range-separated RPA variants for weak intermolecular interactions [J. Toulouse, W. Zhu, 
A. Savin, G. Jansen, and J. G. Angyan, J. Chem. Phys. 135, 084119 (2011)], this works confirms range-separated 
RPAx-S02 as a promising method for general chemical applications. 


I. INTRODUCTION 

There has recently been a revived interest in the random- 
phase approximation (RPA) for the calculation of electron cor¬ 
relation energies of atomic, molecular and solid-state systems 
(see Refs. 1 and 2 for recent reviews). One of the advan¬ 
tages of RPA is its good description of dispersion interactions 
at large separation,^^^ whereas two important disadvantages 
are its poor description of short-range electron correlations® 
and its slow Gaussian basis convergence.These two lim¬ 
itations can be overcome by the range-separation approach 
(see, e.g.. Ref. 10) which allows one to combine a short-range 
density-functional approximation with a long-range RPA-type 
approximation. ^ 

As a result of the increasing interest in RPA, there are now 
several RPA formulations in which RPA equations can be de¬ 
rived, namely the adiabatic connection,^’^^’^^ dielectric ma- 
trix, plasmon ’ formula and ring coupled-cluster dou- 
bles ’ ’ formulations. Moreover, within these formula- 
tions, many variants of RPA {e.g., direct RPA,^^ RPA with ex¬ 
act Hartree-Fock (HF) exchange,RPA with exact Kohn- 
Sham exchange^') can be defined (see also Refs. 16, 17, 32, 
and 33). 

Recently, in a series of papers, we have studied the per¬ 
formance of these RPA variants for closed-shell systems, es¬ 
pecially in the range-separated framework. Within the ring 
coupled-cluster doubles formulation, it was found that a spin- 
restricted closed-shell RPA variant with exchange, first pro¬ 
posed by Szabo and Ostlund (see equation (3.22) of Ref. 28) 
and referred to as Szabo-Ostlund’s second formula (RPAx- 
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S02) in Ref. 16, was one of the best RPA variants for cor¬ 
relation energies, chemical reaction energies, and weak inter¬ 
molecular binding energies for closed-shell systems with or 
without range separation. In view of future applications, 
like the simulation of reactivity in condensed phases, where 
different kinds of intermolecular and intramolecular forces are 
acting simultaneously in a concomitant manner, it is interest¬ 
ing to know whether we can expect a similar degree of reli¬ 
ability in the description of chemical transformations as for 
weak non-covalent interactions. 


The purpose of this work is thus to further test the per¬ 
formance of the RPAx-S02 variant on more general thermo¬ 
chemistry properties, namely atomization energies and reac¬ 
tion barrier heights, now involving open-shell systems. To 
this end, we have extended and implemented the RPAx-S02 
variant, as well as other RPA variants, in a spin-unrestricted 
formalism with and without range-separation. Another im¬ 
portant domain of application of open-shell range-separated 
RPA may consist of the study of transition metal compounds, 
which is out of the scope of the present study and will be the 
subject of forthcoming publications. 


The paper is organized as follows. In Section II, we of¬ 
fer an overview of the RPA variants in the adiabatic connec¬ 
tion and ring coupled-cluster doubles formulations in a spin- 
unrestricted formalism. After giving computational details in 
Section III, we present in Section IV the results of the different 
RPA variants with and without range-separation on datasets of 
atomization energies and reaction barrier heights. A detailed 
statistical analysis of the results is given. Finally, Section V 
contains our conclusions. 
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II. SPIN-UNRESTRICTED RPA CORRELATION ENERGY 
EXPRESSIONS 

We start by giving the RPA correlation energy expressions 
used in this work, using real-valued canonical spin-orbitals in 
a spin-unrestricted formalism for general open-shell systems, 
breaking spatial and 5^-spin symmetry (but not .S^-spin sym¬ 
metry). 

The RPA correlation energy expressions can be derived 
in different/ormu/af/oui (see Ref. 17). One way to appre¬ 
ciate the links between the formulations is to consider the 
adiabatic-connection fluctuation-dissipation theorem expres¬ 
sion for the correlation energy,^^’^® which involves integra¬ 
tions over both a frequency and a coupling constant. An an¬ 
alytical integration over the frequency variable followed by 
a numerical integration over the coupling constant yields the 
adiabatic-connection formulation, while an analytic integra¬ 
tion over the coupling constant followed by a numerical in¬ 
tegration over the frequency yields the dielectric formulation. 
In certain cases, the two integrations can be carried out ana¬ 
lytically which leads to two other formulations: the plasmon- 
formula formulation and the ring coupled-cluster doubles for¬ 
mulation. Depending on the approximations used, the cor¬ 
relation energies derived in these four formulations are in 
general not equivalent. Here we consider only the adiabatic- 
connection and ring coupled-cluster doubles formulations. 

In the adiabatic-connection formulation, the RPA correla¬ 
tion energy is written as the integral over the coupling constant 
a of the trace of the correlation part of a two-particle density 
matrix Pc,a contracted with a matrix Bi of two-electron in¬ 
tegrals at or = 1. Depending whether the calculation of Pc,a 
includes only a direct term (dRPA) or also a HF exchange term 
(RPAx), and whether the two-electron integrals in Bi are non- 
antisymmetrized (I) or antisymmetrized (II), four RPA corre¬ 
lation energy expressions can be defined: the dRPA-I variant 
(also called dRPA or just RPA in the density-functional theory 
literature)^^’^® 


^ da tr [b^ , (1) 

the dRPA-II variant'^ 

£“ = ^‘datr[BfPf/^], (2) 

the RPAx-I variant“d5 

^RPAx-I ^ 1 pRPAxj ^ (3) 

and the RPAx-II variant^^’^^ 


t-RPAx-II 

K 



pRPAxl 
^c,a \ ' 


(4) 


where M„ = (A„ - B^)*/^ (A„ -f B„) (A„ - B,,)'^^ The dRPA 
density matrix is obtained using the following definition 
for the matrices Aq, and for an arbitrary coupling constant 
a. 


jb ^ “ ^ddijSab + a {ib\ai), (6) 

KLj.=«<■') 


where /, j and a, b refer to occupied and virtual spin-orbitals, 
respectively, e, and Ba are the spin-orbital energies, and {ib\a j) 
and {ab\ij) are non-antisymmetrized two-electron integrals. 
Similarly, the RPAx density matrix is obtained using 

the matrices 


(a“).^ ^ “ e,)diAi + {ib\\aj), (8) 

where (//7||ay) = {ib\aj)-{ib\ja} and {ab\\ij) - {ab\ij)-{ab\ji} 
are antisymmetrized two-electron integrals. 

In the ring coupled-cluster doubles formulation, the key 
quantity replacing Pc,a is the matrix of double-excitation am¬ 
plitudes T given by the Riccati equation (at a - 1)^^ 

Bi- hAi T-hTAi h-TBi T = 0 . (10) 

The dRPA-I and RPAx-II correlation energies are then ob¬ 
tained by contracting the amplitudes with two-electron inte¬ 
grals^’ 


^dRPA-I ^ 1 jj. I^gl rpdRPAj ^ (J J) 

£RPAx-II^ Ij^^BlI^RPAxj^ (12) 

where is obtained using the matrices A^ and Bj in 

Eq. (10), and is obtained using the matrices A*/ and 

BJ^ in Eq. (10). We emphasize that the correlation energies 
given by Eqs. (11) and (12) are identical to the ones given 
by Eqs. (1) and (4). Additionally, two other RPA correlation 
energy variants can be defined in this formulation: a dRPA 
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variant with second-order screened exchange (SOSEX) ’ ’ 


c-SOSEX 

Ec 



[b“ , 


(13) 


and a RPAx variant corresponding to Szabo-Ostlund’s second 
formula (S02) (see Ref. 16 and equation (3.22) in Ref. 28) 


£RPAx-S 02 ^ ^RPAxj 


(14) 


In all these cases, Pc,q; is calculated as' 

P,,„ = (A„ - B„)'/2 (A„ - B„)'/2 - 1, (5) 


The expressions of these last two RPA correlation energy vari¬ 
ants are similar but not equivalent to the variants in Eqs. (2) 
and (3). The RPAx-S02 correlation energy expression in 
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Eq. (14) is an obvious extension of the closed-shell Szabo- 
Ostlund’s expression to a spin-unrestricted formalism, and is 
original to this work. We note that there is no obvious exten¬ 
sion of the closed-shell RPAx variant corresponding to Szabo- 
Ostlund’s first formula (SOI) (see Ref. 16 and equation (3.20) 
in Ref. 28) to a spin-unrestricted formalism. At second order 
in the electron-electron interaction, all the RPA correlation en¬ 
ergy variants reduce to the second-order Mpller-Plesset (MP2) 
correlation energy expression,*^ except for the dRPA-I variant 
which reduces to direct MP 2 .^® 

The dRPA-I, dRPA-II, and SOSEX variants are free from 
instabilities by construction. In a spin-restricted closed-shell 
formalism, the RPAx-I and RPAx-S02 variants involve only 
spin-singlet excitations and are thus not subject to triplet insta¬ 
bilities, while the RPAx-II variant involves both spin-singlet 
and spin-triplet excitations and are thus prone to triplet in¬ 
stabilities. Because for spin-unrestricted calculations the S^- 
spin symmetry is imposed (i.e., each spin-orbital has a defi¬ 
nite spin t or j,), the spin-flipped excitations (when the spin- 
orbitals i and a have different spins, or the spin-orbitals j and b 
have different spins) do not contribute to the dRPA-I, dRPA-II, 
RPAx-I, SOSEX, and RPAx-S02 correlation energy variants. 
This can be exploited to reduce the cost of the calculations 
and to avoid that the evaluation of the correlation energy in 
these RPA variants be contaminated by spin-flipped instabil¬ 
ities. On the other hand, the spin-flipped excitations do con¬ 
tribute to the RPAx-II correlation energy variant^** which is 
thus subject to spin-flipped instabilities. Eor this reason, and 
because ofits general poor performance,*® *^’^"* we do not con¬ 
sider the RPAx-II variant in this work. 

RPA calculations can be carried out with and without range 
separation. Eor the full-range case, RPA correlation energies 
are evaluated using orbitals and orbital energies from a HP 
self-consistent calculation. Eor the range-separated case, the 
self-consistent-field (SCP) starting point is a range-separated 
hybrid (RSH) calculation:"*^ 

Ersh = mm {<0|f + + Cl®) + ^hxc [«®]} > (15) 

where <5 is a single-determinant wave function of density n®, 
T is the kinetic energy operator, Vne is the nuclei-electron 
interaction operator, IT^g is a long-range electron-electron 
interaction operator associated with the long-range interac¬ 
tion Wgg(r) = erf (fir) I r, and i® ’■he corresponding jj.- 

dependent short-range Hartree-exchange-correlation density 
functional. The parameter fi controls the range of the sepa¬ 
ration. The RSH energy does not contain long-range corre¬ 
lation, which is added afterwards by evaluating the RPA cor¬ 
relation energies using long-range two-electron integrals and 
RSH orbitals and orbital energies. 

III. COMPUTATIONAL DETAILS 

All calculations have been performed with a development 
version of the MOLPRO 2012 program,"** in which our previous 
spin-restricted implementations of the RPA correlation energy 
expressions have been generalized to a spin-unrestricted for¬ 
malism. We use a straightforward implementation in which 


all the RPA methods scale as where No and Ay are the 

numbers of occupied and virtual orbitals. 

The full-range calculations (post-HP calculations) are la¬ 
belled HP-t-dRPA-I, HP-i-dRPA-II, etc. Eor the range- 
separated calculations, we use the short-range Perdew-Burke- 
Ernzerhof (srPBE) exchange-correlation functional of Ref. 43 
(which is a modified version of the one of Ref. 44), 
and the post-RSH calculations are labelled RSH+dRPA-I, 
RSH-tdRPA-II, etc. The effect of the choice of the short-range 
density functional was studied in a number of previous works 
(see, e.g.. Refs. 14, 18, and 45 ). We do not expect a large de¬ 
pendence of the results on the short-range density functional, 
nevertheless a comprehensive study on the subject would be 
useful. 

We have performed a wide range of RPA calculations on the 
AE 6 and BH 6 datasets"*^ to assess the dependence on the basis 
and on the range-separation parameter yU, as well as the per¬ 
formance of the different methods. The AE 6 dataset is a small 
representative benchmark of six atomization energies consist¬ 
ing of SiH 4 , S 2 , SiO, C 3 H 4 (propyne), C 2 H 2 O 2 (glyoxal), and 
C 4 H 8 (cyclobutane). The BH 6 dataset is a small representa¬ 
tive benchmark of forward and reverse hydrogen transfer bar¬ 
rier heights of three reactions, OH - 1 - CH 4 —> CH 3 + H 2 O, H 
-t- OH —> O -F H 2 , and H - 1 - H 2 S ^ HS + H 2 . We compute 
mean absolute deviations (MADs) as a function of the range- 
separation parameter ^ and as a function of the cardinal num¬ 
ber of Dunning cc-pVXZ basis sets."*® All the calculations for 
the AE 6 and BH 6 datasets were performed at the geometries 
optimized by quadratic configuration interaction with single 
and double excitations with the modified Gaussian-3 basis 
set (QCISD/MG3).®*’ The reference values for the atomiza¬ 
tion energies and barrier heights are the non-relativistic PC- 
CCSD(T)/cc-pVQZ-P12 values of Refs. 46 and 47. 

We have further tested the MP2, dRPA-I, and RPAx-S02 
methods on a set of 49 atomization energies®* referred here 
to as the AE49 dataset (consisting of the G2-1 dataset®^’®® 
stripped of the six molecules containing Li, Be, and Na) and 
on the DBH24/08 dataset®"*’®® of 24 forward and reverse re¬ 
action barrier heights. These calculations were performed 
with the cc-pVQZ basis set, with MP2(full)/6-31G* geome¬ 
tries for the AE49 dataset, and with the aug-cc-pVQZ basis 
set with QCISD/MG3 geometries for the DBH24/08 dataset. 
We have carried out a statistical analysis of the results in the 
form of normal distribution of errors. The reference values 
for the AE49 dataset are the non-relativistic PC-CCSD(T)/cc- 
pVQZ-P12 values of Ref. 56, and the reference values for the 
DBH24/08 dataset are the zero-point exclusive values from 
Ref. 55. 

The calculations labelled dRPA-I, dRPA-II and RPAx-I cor¬ 
respond to Eqs. (l)-(3), where the integral over the coupling 
constant is carried out by a 7-point Gauss-Legendre quadra¬ 
ture. The calculations designated by SOSEX and RPAx-S02 
correspond to Eqs. (13) and (14). The Riccati equations 
are solved iteratively, as described in Ref. 16. Core elec¬ 
trons are kept frozen in all our calculations. Spin-restricted 
calculations are performed for all the closed-shell systems, 
and spin-unrestricted calculations for all the open-shell sys¬ 
tems. Eor a few systems, the full-range RPAx-I integrand in 
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FIG. 1. Mean absolute deviations for the AE6 and BH6 datasets 
as functions of the range-separation parameter yu for range-separated 
calculations using the cc-pVQZ basis set, the srPBE functional and 
different post-RSH long-range correlation methods (MP2, dRPA-I, 
dRPA-II, RPAx-I, SOSEX, and RPAx-S02). The case /z = oo corre¬ 
sponds to full-range post-HE calculations. The reference values are 
the non-relativistic FC-CCSD(T)/cc-pVQZ-F12 values of Refs. 46 
and 47. 


FIG. 2. Mean absolute deviations for the AE6 and BH6 datasets as 
functions of the cardinal number X of the Dunning cc-pVXZ basis 
sets for full-range post-FlF and range-separated post-RSH (srPBE 
functional and ji = 0.5) calculations with different correlation meth¬ 
ods (MP2, dRPA-I, dRPA-II, RPAx-I, SOSEX, and RPAx-S02). The 
reference values are the non-relativistic EC-CCSD(T)/cc-pVQZ-E12 
values of Refs. 46 and 47. 


Eq. (3) diverges in the vicinity of tr = but in practice 
in our calculations since the quadrature abscissae do not ex¬ 
ceed a — 0.975 this divergence is avoided. A more detailed 
study of instabilities will be shown elsewhere. In practice, we 
did not encounter instabilities for the full-range HF-hRPAx- 
S02 method, provided of course that the true SCF minimum 
was found. Furthermore, we did not observed any instabilities 
for the range-separated RSH-i-RPAx-I and RSH-i-RPAx-S02 
methods. 


IV. RESULTS AND DISCUSSION 

Throughout the paper, the following color and point symbol 
code is used: MP2 results are shown in black and their symbol 
is a circle; dRPA-I figures are displayed in blue, their symbol 
is a square; and RPAx-S02 data are drawn in red, their symbol 
being a triangle. When a distinction needs to be done between 
full-range calculations and range-separated calculations, the 
former are empty symbols and the latter are solid symbols. 


Figure 1 shows MADs for the AE6 and BH6 datasets as a 
function of the range-separation parameter fi, from jj. - 0.3 
bohr *, where all methods nearly coincide since the long- 
range correlation contribution is small, to full-range calcu¬ 
lations (yu = oo), where the methods yield very different re¬ 
sults. Range separation greatly reduces the MADs on the two 
datasets for all the methods. For the AE6 dataset, the usually 
used value of yu = 0.5^’ for the range-separation parameter 
yields results reasonably close to optimal for all the methods, 
with MADs around 5 kcal/mol. In the case of the BH6 dataset, 
although the value of ji — 0.6 gives overall the best results, 
the MADs obtained with /r = 0.5 lie within 1 kcal/mol of the 
minimal MADs. Consequently, we choose to use the value of 
yu = 0.5 in all the following. Note that for this value of yu all 
the range-separated methods give about the same MADs. 

Figure 2 shows the MADs for the AE6 and BH6 datasets 
as a function of the cardinal number X of the Dunning cc- 
pVXZ basis sets (X=2,3,4). For the AE6 dataset, the overall 
gain in MAD when going from cc-pVDZ to cc-pVQZ calcu¬ 
lations is around 25 kcal/mol for all full-range methods except 
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FIG. 3. Distribution of the errors (in kcal/mol) obtained with full- 
range (post-HF) and range-separated (post-RSH) calculations on the 
AF49 dataset using the cc-pVQZ basis set and yu = 0.5. The bins are 
the distributions of the actual errors and the curves are fitted Gaus¬ 
sian distributions. The reference values are the non-relativistic FC- 
CCSD(T)/cc-pVQZ-F12 values of Ref. 56. 


FIG. 4. Distribution of the errors (in kcal/mol) obtained with full- 
range (post-HF) and range-separated (post-RSH) calculations on the 
DBH24/08 dataset using the aug-cc-pVQZ basis set and fi = 0.5. 
The bins are the distributions of the actual errors and the curves are 
fitted Gaussian distributions. The reference values are taken from 
Ref. 55. 


TABLE I. Mean absolute deviations (in kcal/mol) for the AE49 
and DBH24/08 datasets for range-separated calculations using the 
srPBE functional and several values of the range-separated param¬ 
eter (41 = 0.3, /J = 0.5 and = 0.7) with different post-RSH long- 
range correlation methods (MP2, dRPA-I, and RPAx-S02). The case 
yu = 00 corresponds to full-range post-HF calculations. The basis 
sets are cc-pVQZ for the AE49 dataset and aug-cc-pVQZ for the 
DBH24/08 dataset. The reference values are the non-relativistic FC- 
CCSD(T)/cc-pVQZ-F12 values of Ref. 56 for the AE49 dataset and 
the values of Ref. 55 for the DBH24/08 dataset. 


AE49 DBH24/08 


= 

0.3 

0.5 

0.7 

00 

0.3 

0.5 

0.7 

00 

MP2 

6.34 

5.09 

7.00 

5.63 

8.28 

2.94 

3.45 

6.17 

dRPA-I 

6.38 

6.80 

12.09 

23.23 

8.29 

3.01 

3.89 

6.95 

RPAX-S02 

6.50 

4.06 

5.08 

12.20 

5.12 

2.83 

3.27 

6.40 


for RPAx-S02 whose MAD increases (the mean deviation 
changes sign, not shown), while the gain for range-separated 
calculations is around 10 kcal/mol. Similarly, the gain when 
going from cc-pVTZ to cc-pVQZ is about 5 kcal/mol for full- 
range calculations and is negligible for range-separated calcu¬ 
lations. The full-range calculations are thus not yet converged 


with the cc-pVQZ basis set while the range-separated calcu¬ 
lations can be considered as converged with the same basis 
set. This clearly demonstrates that range-separated methods 
are much less basis-set dependent than full-range methods for 
calculations of atomization energies. These observations were 
expected given the demonstrated exponential convergence of 
long-range correlation energies with respect to the cardinal 
number.^^ For reaction barrier heights, both the full-range and 
range-separated methods are relatively weakly dependent on 
the basis set. 

We now discuss results for the full-range and range- 
separated MP2, dRPA-I, and RPAx-S02 methods on the 
larger AE49 and DBH24/08 datasets. Table I shows MADs 
for selected values of the range-separation parameter p. For 
the full-range case - 00 ), the MADs that we obtain on the 
AE49 dataset for HF+MP2 and HF+dRPA-I (5.63 kcal/mol 
and 23.23 kcal/mol, respectively) are in good agreement with 
the MADs on the G2-I dataset reported with the same meth¬ 
ods with a different basis set in Ref. 59 (5.9 kcal/mol and 21.7 
kcal/mol, respectively). 

For range-separated calculations, the results confirm that, 
among the values tested, p = 0.5 yields the smallest MADs 
for all the methods for the two datasets. For this value 
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FIG. 5. Representation of the means m and standard deviations cr (in kcal/mol) of the Gaussian distributions fitted to the distribution of errors 
seen in Figures 3 and 4 for full-range (post-HF) and range-separated (post-RSH) calculations on the AE49 and DBH24/08 datasets. 


of /i, the RSH-t-RPAx-S02 method gives the smaller MAD 
of 4.06 kcal/mol for the AE49 dataset, slightly better than 
RSHh-MP 2 which gives a MAD of 5.09 kcal/mol. For the 
DBH24/08 dataset, at = 0.5, the RSH+MP2, RSH-t-dRPA- 
I, and RSH-i-RPAx-S02 methods all give similar MADs of 
about 3 kcal/mol. These results confirm the representative¬ 
ness of the AE6 and BH6 datasets and justifies the subsequent 
analysis of the range-separated MP2, dRPA-I, and RPAx-S02 
methods at yU = 0.5. 

We note that full-range RPA calculations are often per¬ 
formed as post-Kohn-Sham calculations, i.e. using Kohn- 
Sham orbitals and orbital energies. We choose instead to per¬ 
form full-range RPA calculations as post-HF calculations be¬ 
cause it corresponds to the limit of RSH-i-RPA calculations 
for fj. —> oo. It turns out that the MAD that we obtain with 
HF-i-RPAx-S02 on AE49 (12.20 kcal/mol) is similar to the 
MADs reported for post-PBE dRPA-I calculations on the G2- 
I dataset (13.3 kcal/mol^® or 10.2 kcal/mol®° depending on the 
basis sets). 

The detailed results of the calculations on the AE49 and 
DBH24/08 datasets for the full-range and range-separated 
MP2, dRPA-I, and RPAx-S02 methods at/j - 0.5 are given in 
the supplementary material.®* These data are analyzed in Fig¬ 
ures 3 and 4 using distributions of errors. The bins correspond 
to the actual distributions of the errors and the curves are fit¬ 
ted Gaussian distributions of mean m and standard deviation 
cr. For the AE49 dataset, full-range MP2 gives a small mean 
error of m = -0.57 kcal/mol but with a standard deviation of 
cr = 7.02 kcal/mol, while full-range dRPA-I strongly underes¬ 
timates (with a mean of m = -23.23 kcal/mol) and full-range 
RPAx-S02 overestimates (m = 9.33 kcal/mol) the atomiza¬ 
tion energies. When going from full-range to range-separated 
calculations, the three methods give much narrower distribu¬ 
tions of the errors, and for the case of dRPA-I and RPAx-S02 
much smaller mean errors of m = -6.43 and m - -2.90 
kcal/mol, respectively. 

For the DBH24/08 dataset, the three full-range methods 
tend to overestimate the reaction barrier heights. Range sepa¬ 
ration does not change much the standard deviations (which 


range from 4.81 kcal/mol to 7.08 kcal/mol), but it greatly 
reduces the mean errors which are m = 0.20, m = 0.55, 
and m = -0.19 kcal/mol for RSH-t-MP2, RSH-t-dRPA-I, and 
RSH-i-RPAx-S02, respectively. 

Finally, the means m and standard deviations cr of the fitted 
Gaussian distributions are reported in Figure 5, which gives a 
way of easily assessing the performance of the full-range and 
range-separated methods by their position in the (m, cr) plane. 
If we define the best method as the one closest to the point 
(0,0) in this plane, the best method is RSH-t-RPAx-S02 for 
atomization energies while RSH-)-MP2 and RSH-i-RPAx-S02 
are about equally good for reaction barrier heights. 

Given that for atomization energies and reaction barrier 
heights RSH-t-MP2 is about as accurate as RSH-t-RPAx- 
S02 while being less expensive, one could conclude that 
RSH-)-MP2 is a better method. However, RSH4-RPAx-S02 
is significantly more accurate that RSH-I-MP2 for weak inter- 
molecular interactions (see Ref. [ 16 ]), hence RSH-t-RPAx- 
S02 appears as a more systematic method for various types 
of applications. 


V. CONCLUSION 

We have implemented several RPA correlation energy vari¬ 
ants in a spin-unrestricted formalism with or without range 
separation, and tested them on thermochemistry datasets of 
atomization energies and reaction barrier heights. Range- 
separation greatly improves the accuracy of all the RPA vari¬ 
ants. Specifically, the RSH-t-RPAx-S02 variant is among the 
most accurate range-separated RPA variant with mean ab¬ 
solute deviations of about 4 kcal/mol for atomization ener¬ 
gies and about 3 kcal/mol on reaction barrier heights. Since 
RSH-(-RPAx-S02 had already been shown to be among the 
most accurate range-separated RPA variant for weak inter- 
molecular interactions,*® this work confirms RSH-)-RPAx- 
S02 as a promising method for general chemical applications. 
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